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Abstract 

Various problems concerning the geometry of the space u*{Ti.) of Hermitian operators on a 
Hilbert space Ti. are addressed. In particular, we study the canonical Poisson and Riemann- Jordan 
tensors and the corresponding foliations into Kahler submanifolds. It is also shown that the space 
ViTi) of density states on an n-dimensional Hilbert space Ti. is naturally a manifold stratified 
space with the stratification induced by the the rank of the state. Thus the space 7?'^(7i) of rank-fc 
states, fc = 1, . . . , n, is a smooth manifold of (real) dimension 2nk — — 1 and this stratification 
is maximal in the sense that every smooth curve in T>{Ti), viewed as a subset of the dual M*(7i) to 
the Lie algebra of the unitary group UiTi), at every point must be tangent to the strata V'^iTi) it 
crosses. For a quantum composite system, i.e. for a Hilbert space decomposition TL — Ti} ® Ti^, 
an abstract criterion of entanglement is proved. 

1 Introduction 

Dirac's approach to Quantum Mechanics uses a Hilbert space as a fundamental object to start with, 
motivating the linear structure with the superposition principle necessary to describe phenomena like 
those of interference J]. Born's probabilistic interpretation requires the use of a Hermitian inner 
product to deal with normalized states, therefore the physical identification of states in the Hilbert 
space leads to the requirement that (pure) states of a quantum mechanical system are described by 
elements of the complex projective space (one-dimensional subspaces of a separable complex Hilbert 
space Ti.). By means of the Hermitian structure on Ti it is possible to define a binary product on 
the pure states PTi El E] • The physical interpretation of this binary operation is given in terms 
of probability transition from one state to another. On this space PTi, bijective maps which preserve 
the transition probability are necessarily projection of unitary or anti-unitary transformations on the 
original Hilbert space, this statement is the main content of Wigner's theorem 
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More likely, due to this "equivalence" between the two descriptions (on TL and on P7i), physicists 
have barely paid any attention to the geometrization of Quantum Mechanics, i.e. to introduce a 
"tensorial description" in such a way that non-linear coordinate transformations could be performed, 
notably exception obviously do exist and we provide a partial list of references ^21 • The recent great 
interest in the foundational aspects of Quantum Mechanics motivated by the use of entanglement as 
a resource for quantum information and quantum computing has boosted a more deep study of many 
fundamental aspects, for instance the possibility to have a binary composition of pure states without 
the use of the Hilbert space linear structure [51 the possibility to have a non- linear Quantum 
Mechanics ^JEI, more generally the possibility to have non-linear transformations among states. 

The possibility of non-linear transformations may turn out to be quite useful in the classification 
problem of separability and entanglement because these properties are not preserved by taking linear 
combinations. Moreover, an appropriate description of atomic phenomena involving polarization, spin 
orientation and angular correlations, requires that we go beyond pure states in the description of 
quantum systems. This larger family of states was introduced by von Neumann as dual objects with 
respect to the quantum observables, they constitute the set of density states and an early, physically 
motivated, review was written by U. Fano 

Again, for these states a proper mathematical setting is provided by the dual space of the Lie algebra 
of the observables, with respect to the coadjoint action of the unitary group. Density states emerge 
as elements of the coadjoint orbits passing trough some special elements in the dual of the Cartan 
subalgebra. The mathematical context of coadjoint orbits is quite well known to those physicists 
involved with geometric quantization and the field was widely studied in the seventies by Kostant, 
Kirillov, and Souriau [T^. 

Each coadjoint orbit bears a natural differential structure. Observe, however, that the spectrum of 
the state does not change along the orbit of the unitary action. From the point of view of quantum 
evolution it corresponds to the situation of an isolated system, when all interactions with the envi- 
ronment are negligible, so there is no dissipation and the evolution is unitary. In many cases this is 
only a very exceptional situation, very rarely adequately corresponding to the physical reality. On the 
other hand, it is a priori not clear that the whole set of density states, i.e. a union of coadjoint orbits 
of the unitary action of different dimensionality, possesses a natural differential structure. Exhibiting 
such a structure in terms of local coordinates and/or via a general geometric construction of a smooth 
stratification of density states is thus of great interest when investigating dissipative systems. 

Density states form a convex subset of the set of Hermitian operators on Ti. Some properties of 
these convex body attracted recently an attention ^H]. It is thus legitimate to ask about "the shape" 
of the set of density matrices, in particular about the smoothness properties of its boundary. In the 
simplest case of the two-dimensional Ti, the density matrices form the three-dimensional unit ball with 
a smooth boundary - the two-dimensional unit sphere comprising all pure states. But this situation 
is exceptional - in higher dimensions the boundary does not consists exclusively of pure states, it is in 
addition not smooth. 

The space of density states carries additional structures with respect to those available on the space 
of coadjoint orbits of general Lie groups because they are related to the unitary group and therefore 
additional structures are available. Moreover the need to consider composite quantum systems, tensor 
products of the spaces associated with a choice of subsystems making up the whole system, will bring 
up novel problems which will require further investigations. 

All these various considerations have convinced us that a review of these mathematical aspects along 
with the identification of the novel emerging problems may be useful to those people interested in the 
application of quantum mechanics to quantum information and are not at home with the geometrical 
background required. A recent book by Chruscihski and Jamiolkowski deals with geometrical 
aspects of quantum mechanics, these authors however are primarily concerned with the application 
of these methods to describe the geometric phase ^7]. At this point one should also point to the 
paper '18' in which, in connection with geometric phase and parallel transport along mixed states, the 
geometry of the manifold of density matrices as a stratified space, was discussed along slightly different 
lines than in the present paper (cf. Section 3 below). 

In the present paper the Hilbert space Ti. will be assumed to be of finite dimension n in order to 
make the differential geometry expressible in local coordinates classical. The reader will understand 
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that passing to an infinite-dimensional H. (i.e. a differential geometry of a Banach (or a Hilbert) 
manifold) is straightforward, to this aim we will try to use coordinate-free expressions, which serve in 
both cases, as much as possible. The paper is organized as follows : 

In section 2 we start with presenting the Kahler structure on the Hilbert projective space PTL of 
pure states obtained from the standard Hermitian product on Ti via the momentum map associated 
with the Hamiltonian action of the group U{Ti.) of unitary transformations of Ti.. In this picture the 
pure states form just an orbit in the dual space u*{H) of the unitary Lie algebra u{Ti) of the group 
UiTi). Because of the nondegeneracy of the canonical invariant scalar product on uiTi.) we have a 
canonical identification of u*{'H) with u{H) which makes the geometry of u*{H) very rich. We decided 
to interpret u* {Ti) as the space of Hermitian operators on H which makes possible to understand the 
density states as a subset of u*{H). 

Consequently, in sections 3 and 4 we present the density states as a convex body ViTi) in u*{TL) 
which is a family of some C/(7i)-orbits and, as we will show later, also orbits of a particular action 
of the group GL{Ti) of invertible complex linear operators on TL. We show that I?(7i) is naturally a 
manifold stratified space with the stratification induced by the the rank of the state. Thus the space 
I''^(7i) of rank-fc states. A: = 1, . . . , n, is a smooth manifold of (real) dimension 2nk — ~ 1 and this 
stratification is maximal in the sense that every smooth curve in ViTi.), viewed as a subset of the dual 
u*{'H) to the Lie algebra of the unitary group U(H), at every point must be tangent to the strata 
'D'^(H) it crosses. 

Section 5 is devoted to the geometry of u* {Ti.) , to a global description of the Kahlerian structure 
of U (7i)-orbits by means of the canonical Poisson and Riemann-Jordan tensors. These Kahlerian 
structure are well-known in algebraic geometry and can be easily generalized to analogous structures 
on general flag manifolds. The point which should be stressed here is that the geometry we develop is 
canonical, that it does not depend on the matrix form of an operator and the C/(7i)-orbits are treated 
as a collection rather than each orbit separately. 

In the last section we investigate a Hilbert space decomposition Ti = Ti} (g) TL^ which is usually un- 
derstood as corresponding to a quantum composite system. We present an introduction to the problems 
of separability and entanglement together with an abstract scheme for measurement of entanglement. 

Geometry of composite quantum systems was investigated in the literature from several points of 
view. First, it is of importance to distinguish classes of states which are equivalent under a restricted 
set of unitary transformations (dubbed local transformations in the physical literature), namely those 
which belong to the same orbit of U{Ti}) x U{Ti?). From the physical point of view all states on the 
same orbit contain an equal amount of quantum correlations between the subsystems, i.e., these can 
not be influenced by operations performed separately on each subsystem. 

In order to characterize uniquely an orbit (i.e. a class of locally equivalent states) one can try to find 
a complete set of U {Ti^) x U (7i^)-invariant functions on T>{H), such that the values of all functions at 
p £ ViTi) characterize uniquely the orbit through p |19ll2()j . The task can be effectively completed only 
for low-dimensional systems - in fact, only in the the case dim7i^=dim7i^ — 2 the explicit results were 
found ,21^. The same is true for multipartite composite systems i.e. when Ti, = TC^ (gi Ti.^ (g • ■ • (g Tl^ . 
Here also the explicit results are known for X up to 4 and dimW ~ 2, i ~ 1, . . . ,4 \22i I23j . 

Other (partial) characterization of local orbits is provided by their dimensions. These were inves- 
tigated in 1^24, ,2^ and in j22j all orbits of submaximal dimensionality in the case dim?i^=dim?i^ = 2 
were explicitly identified and enumerated. The similar task of finding dimension of the local orbit 
through an arbitrary p in the case of higher dimensional systems was never achieved. A much modest 
goal of determining dimensions and topology of local orbits stratifying the set of rank one (pure) states 
'D^{Ti} ® Ti^) was, however, completed for arbitrary finite-dimensional Ti} and Ti^ 27 . 

The sets of pure states in two- and three-partite systems with dimTi* = 2 can be identified with, 
respectively, unit seven- and fifteen- dimensional spheres S*" and S^^. In both cases there exist the 
Hopf fibrations S'' — > and S^^ — > which were used to investigate the geometry of pure states in 
|28l f29" '^D', TT , whereas multipartite pure states were treated in using Segre variety. 

Although in the present paper we limit ourselves to investigation of two-partite composite system, 
we would like to point out recent achievements in geometric characterization of entangled pure states 
of multipartite systems. When investigating entanglement in multicomponent system one aims at 
discriminating among different classes of entanglement, defined as different equivalence classes under 
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appropriate group of transformations preserving entanglement properties. The goal can be achieved by 
identification of the so called entanglement monotones, i.e. measures of entanglement which are invari- 
ant under considered transformations. Construction of such invariants based on Pliicker coordinates 
on Grassmannians, naturally appearing when considering pure states of multicomponent systems, were 
presented in [IVA\ and j84j . The geometry of three-qubit pure states entanglement was recently inves- 
tigated in [35|. where geometric description of different classes of entanglement was given in terms 
of submanifolds of the so-called Klein quadric - a special quadric embedded in the five-dimensional 
complex projective space. 

Of special interest is the set of separable states (defined in Sec. 6), as those which, from the physical 
point of view, do not carry any quantum correlations. From the construction they form a convex subset 
in T>{H^ ®Ti?). Only in the case of dim7Y^=dimH^ = 2 and dmiH} = 2 and dimTC^ = 3 (or vice versa) 
there exist effective criteria which allow to discriminate separable and nonseparable (entangled) states. 
As a consequence only in these low-dimensional case one can relative easily investigate the geometry 
of the boundary of the set of separable states |37j . 

2 Kahler structure on the Hilbert projective space 

Let Ti be an n-dimensional Hilbert space with the Hermitian product (x, y)-H being, by convention, 
C-linear with respect to y and anti-linear with respect to x. The unitary group U{7i) acts on TC 
preserving the Hermitian product and it consists of those complex linear operators A e gKji) on H 
which satisfy AA^ = I, where A^ is the Hermitian conjugate of A, i.e. 

{Ax,y)H = {x,A''y)-H- 

The geometric approach to Quantum Mechanics is based on considering the realification T^r of Ti as 
a Kahler manifold {Tin, J, g,Lu) with canonical structures: a complex structure J : THr —^ TT^r, a 
Riemannian metric g, and a symplectic form uj. The latter come from the real and the imaginary parts 
of the Hermitian product, respectively, g = 5R((-, ■)?{): ^ — ^{i', ■)h)- After the obvious identification 
of the tangent bundle TTYr with T^r x Hr, all these structures are constant structures induced from 
H: 

J{x) = i-x, g{x, y) + i- uj{x, y) = (x, 

We have obvious identities 

= -/, u;{x,Jy)=g{x,y), g{Jx, Jy) ^ g{x,y), u;{Jx, Jy) ^ uj{x,y). 

The tensors g and to being non-degenerate have their inverses: the contravariant metric tensor G — g^^ 
and the Poisson tensor = oj^^ . They form together a Hermitian product 

on the dual real Hilbert space TC^ equipped with the dual complex structure J* . Using the identification 
of with TYr via the metric tensor the latter can be interpreted as a contravariant complex tensor 
on Wr. This tensor induces two real brackets of smooth functions on Wr: {/, /i}g = G'(d/, d/i) and 
{f,h}^ = n(d/, d/i). The first one is the 'Riemann- Jordan' bracket associated with the contravariant 
version of the metric tensor g and the second is just the symplectic Poisson bracket associated with w. 
Of course both brackets can be extended to complex functions by complex linearity and give rise to 
the 'total' bracket 

{/, h}H = (d/, dh)n' = {/, h}g + i- {/, h}^. (1) 

Fixing an orthonormal basis (cfe) of TL allows us to identify the Hermitian product {x,y)-H on TL 
with the canonical Hermitian product on C" 

(a, 6)c" = Ofefofc (2) 

(we use the convention of summation on repeated indices), the group U{Ti.) of unitary transformations 
of "H with U{n), its Lie algebra u{'H) with u(n), etc. In this picture {ajkY — (a¥7) a-nd {T^T)jk = 
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(aj, afe), where afe = (tjk) € C" are columns of the matrix T = [tjk)- The choice of the basis induces 
(global) coordinates {qk,Pk), k = 1, . . . , n, on Wr, by 

{ek,x)n = {qk + i-Pk){x) 
in which dq^ is represented by efe and dp^ by i ■ e^. Hence the complex structure reads 

J = dpk® dg/c - dq^ ® dpfe, 

the Riemannian tensor 

9 = {dqk (8> d^fe + dpk (8> dpk) = ^{dqk V dgk + dpk V dpk) 

and the symplectic form 

u) = dqk A dpk, 

where x\/y = xC^y + y^xis the symmetric, and xAy = x^y — y0xis the wedge product. In 
complex coordinates Zk = qk + i ■ Pk one can write the Hermitian product as the complex tensor 

{■,-)n = dzk^dzk- 
The contravariant tensor G + i ■ Ct has the form 

G + i-n = {dq^(g>dq^+ dp^ (g) 9pJ + i ■ (9,, (g) dp^ - dp^ ® dq^) 
or, in complex coordinates, 

G + i-n = {dq^-i-dp^)^{dq^+i-dp^)= Ad^^ O 

In other words, 



^dj_dl^^dj_dh_^ 
^ dqk dqk dpk dpk ' 

dqk dpk dpk dqk ' 

and 

Ifh} =4^ — 
" dzk dzk ' 

Every complex linear operator A e gl{H) on H induces the quadratic function 

fA{x) = ^{x,Ax)n- 
The function is real if and only if A is Hermitian, A = A^ . 

One important convention we want to introduce is that we will identify the space of Hermitian operators 
A = A^ with the dual u*{Tl) of the (real) Lie algebra w(H), according to the pairing between Hermitian 
A e u*{7i) and anti- Hermitian T e u{H) operators 

{A,T)='-.Tr{AT). 

The multiplication by i establishes further a vector space isomorphism u{H) B T iT € u*{H) 
which identifies the adjoint and the coadjoint action of the group J7(7Y), Ad;7(T) = UTU^ . Under this 
isomorphism u*{H) becomes a Lie algebra with the Lie bracket [A,B] = j[A,B]-, where [A, = 
AB — BA is the commutator bracket, equipped additionally with the scalar product 

{A,BU,=^Tr{AB) (3) 



5 



and an additional algebraic operation, the Jordan product [A, = AB + BA. The scalar product is 
invariant with respect to both: the Lie bracket and the Jordan product (or bracket) 

= (4) 

+ = + (5) 

and it identifies once more u*{TC) with its dual, 

u*(n) 3 A^ A=-Ae u(H), 
so vectors with covectors. Under this identification the metric Q correspond to the invariant metric 

(l,B}„ = iTr(AS) (6) 
on u{l-C) which can be viewed also as a contravariant metric on u*{'H). 

For a (real) smooth function / on let us denote by gradf and Hamf the gradient and the Hamil- 
tonian vector field associated with / and the Riemannian and the symplectic tensor, respectively. In 
other words, g{-, gradf) = df and lu{-, Harrif) = df or gradf = G(d/, •) and Haruf = Q{df, •). Note 
that any A e gl{Ti-) induces a linear vector field A on TL by A{x) = Ax. 

Lemma 1 For Hermitian A we have 

gradf ^ — A and Hamf^ — i A. 
Proof. If (•, •) denotes the pairing between vectors and covectors then 

(d/A(x), y) = ^{{y, Ax)n + {x, Ay}n) = 5R((2/, Ax)h) 
— g{a, Ax) — uj{y,iAx). 



Corollary 1 For all A,B€ glCH) we have 

{/a, fsju = f2AB- (7) 

In particular, 

{fAifsjg = Jab+ba, (8) 

{/ylj/sluj = f-i(AB-BA)- (9) 

Proof. For Hermitian A, B we have 

{fA,fB}nix) = gigradAix),gradB{x)) +i ■ uj{Hamf{x),Hamg{x)) 

= g{Ax, Bx) + i ■ uj{iAx, iBx) = {Ax, Bx)u = {x, ABx)u — '^fAB{x). 

But 2AB ^ {AB + BA) + i{-i{AB~BA)), where AB + BA = [A,B]+ and -i{AB-BA) = -i[A,B]^ 
are Hermitian, thus f[A.B]+ ^-nd ^-re real, so the thesis holds for Hermitian A, B. For general 

A, B it follows by complex linearity. I 

The unitary action of U{7i) on Ti is in particular Hamiltonian and induces a momentum map 
M • ''^R ^ u*{7i). The fundamental vector field associated with iA G u{H), where A G u*{TL) is 
Hermitian, reads lA, since 

|t=o exp(--A)(a;) =iA{x). 
dt I 
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The Hamiltonian of iA is Ja, so the momentum map is defined by 

{fi{x), -A) = fA{x) = l-{x,Ax)-H- 
I 2 

But by our convention 

(^(.t), \a) = '-T,{^,{x)\a) = \IT{^^{x)A), 

so that Tr(/i(x)yl) — {x, Ax)-h and finally, in the Dirac notation, 

f,{x)^\x){x\. (10) 

Note that for A being Hermitian Ja is the puUback /a = = ^ ° where A — {A, •)„• = 

£ u{'H). The linear functions A generate T*u*(7i), so that © and mean that the momentum 
map /i relates contravariant tensors G and 51 on respectively, with the linear contravariant tensors 
R and A on u*{Ti.) corresponding to the Jordan and Lie bracket, respectively. The Riemann- Jordan 
tensor i?, defined in the obvious way, 

= + = ^TriaAB + BA)), (11) 

is symmetric and the tensor 

A(0(1, B) = [A, = ^Tr(e(Ai? - i?A)), (12) 

is the canonical Kostant-Kirillov-Souriau Poisson tensor on u*{Ti,). They form together the complex 
tensor 

(i? + ^ • A)(e)(A, B) = 2(e, = Tii^AB) (13) 

and the momentum map relates this tensor with the dual Hermitian product: 

^4G + i-n) ^ R + i- A. (14) 

Example. For TL ~ consider an orthonormal basis in u*{2) consisting of 

"-{\ \)- ^-{l ^-{\ I)' o) 

and the associated coordinates u, x, y, z, where u{A) — iTr(UA), etc. In these coordinates the Poisson 
tensor reads 

A = 2{zd^ Ady + xdy A 9^ + yd;, A d^) 
and the Riemann- Jordan tensor reads 

i? = 5„ V {2xdj, + 2ydy + 2x8^) + u{du dy, + V + dy V dy + V d^). 

The rank of A('u, x, y, z) is if + + = and 2 if + + > 0. The rank of R{u, x, y, z) is 
at (m, X, y, z) = 0, it is 2 for u = and a;^ + + > 0, it is 3 for + + z^ = > 0, and it is 4 
for + j/^ + z^ ^ > 0. 

The image /x('H \ {0}) is the cone 

V^{n) = {\x){x\: x^O} 

of non- negatively defined Hermitian operators ^ =| a;)(a;| of rank 1. The operator ^ is proportional 
to the 1-dimensional projection so = where ||^|| — ||a;|p is the operator norm of ^. 

The manifold V^{H) is foliated by C/(H)-coadjoint orbits being complex projective spaces 'DKH) = 
{\x){x\: \\x\\ = r}, r > 0. In particular, the momentum map image of the (2n-l)-dimcnsional sphere 



7 



Sn = {x & H : ||a;p = {x,x)-h = 1} is the complex projective space T>^{H) = {\ x){x\: \\x\\ = 1} of 
1-dimcnsional projectors. 

The coadjoint orbits O in u*{J-L) possess canonical symplectic forms r]''-' which build together the 
Poisson structure A. These forms, as the inverses of Ajo, are characterized by 

77f([AC],[i3,?]) = = -(^,[AB])«-- (15) 

Indeed, the vectors [A.X] = form the tangent space of the ?7(H)-orbit through O and rf is 

the inverse of A^o- Due to invariancc of the scalar product on u*(Ti.): 

A^AB) = = {[^,A],B)^,^u) = m,A],B). (16) 

Hence #Aj(A) = [^,A] and 

#vmA]) = {#A^)-\[^,A])=A, 



so 

VI{[A,^],[B,^]) = {#vm^mB,^])=-{A,[B,^]) 

= -(A,[B,e])„.(„) = -(C,[A,B])„.(„). 

The image #i?(T*u*(H)) of the tensor R is not an involutive (generalized)distribution, so its inverse 
a = can be understood only as a 'partial' covariant tensor on u*(TC), i.e. as a 'partial symmetric 
2-form' which at ^ G u*{H) is defined only on vectors from #R^{T'^u* {H)). There is a completely 
analogous characterization of the tensor a to that of r]. Both characterizations we can summarize as 
follows. 

Proposition 1 (a) The symplectic form rP on the U{H)-orbit O is characterized by 

nfi[A,^],[B,^]) = = -{^,[A,B]U., (17) 

where A,B £ u* (H) are arbitrary Hermitian operators. 

(h) The 'partial tensor' a on u*(Ti.) is characterized by 

a^i[A,^]+,[B,C] + ) = = {^,[AB]+)u'. (18) 

where A,B £ u* {H) are arbitrary Hermitian operators. 

Let us observe that the 'partial tensor' a, when restricted to any I>^(7i), induces a Ricmannian 
structure which, together with the symplectic structure r]^ = rj^^^^\ induces a Kahler structure. 

Proposition 2 (a) The tensor cr'' being the restriction of the partial tensor a to the U(Ti.)-orbit 
VKH) through ^ = /x(a;), r^ = ||^||, is proportional to the original scalar product on u*{'H): 

al{[A,^],[B,^]) = (19) 

(b) The {l,l)-tensor J on P^{H), J^{A) = satisfies = — J and induces a complex 

structure "^J on every V], (Ti.) . Moreover, 

r,l{[A,iVJd[Bm=<ym^^UB,i]), (20) 

and 

riim[AA]VJd\B.m)=vl{[AmB,i]), (21) 
i.e. {V].{T-C),^ J ,cr^ ,r]^) is a Kahler manifold for each r > 0. 
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Proof. Observe first that, due to the Leibniz rule, 



Then, in view of lfTH|l , 



But 



= tt{[a, e] o [B, e] - [A, e] o e o [B, e] + e ° [b, o [a, c]) 

= Tr([Ae']o[S,e]) = ||e||Tr([Ae]o[B,e]) 

that proves ifT^ . 

To prove that is a. complex structure on every orbit, let us recall that = IICIIC- Passing 
to ^' = 6/11^11 if necessary, we can assume for all the further calculations that = 1 so that 
J^iA) = [A,^]. Hence, 

[p,e],e],e] - ~-ue-2^Ac+eA),^]- = -^[{Ae-eA] = -[a,^] (22) 

and (|19|l follows. Moreover, since vectors [A,£] form the tangent space T^VKn), ^ shows that 
J reduced to 2?^(7i) is an almost-complex structure ""JZ. We shall show that the Nijenhuis torsion 
of ^Jj vanishes, so the structure is integrable. To do this, we must show that the distribution in 
the complexified tangent bundle TX'^(7Y) C which corresponds to eigenvectors of complexified ^J' 
with the eigenvalue i is involutive. But this distribution is generated by complex vector fields T for 
T e gliU), where T{^) = ^T(l - 0- Indeed, 

Jd^Tii - 0) = - 0, = ]{^T{i - - eT{i - 0) = * • ^T{1 - 

and this is a generating set due to the decomposition 

T = im + (1 - 0^(1 - 0) + (1 - + - 

into eigenvectors of J with eigenvalues 0, —i, and i, respectively. The bracket of vector fields [Ti, T2]vf 
reads 

+ eTiCTi(i - = am - T2ri)(i - = m,T2]-) 

that proves involutivity. 

Finally, it is sufhcient to combine ifTTjl and (fH^ to get f^ . Then 

77^^C-j-^([A,0),'-J«([i?,0)) = a^^C-Jc([A,0),[i3,0) = ('■^«([AO),[s,0)n*(«) 

= (P,0,e],[s,0)„*(«) = -([P,e],^],?],i?)«-(H) 

= (P,0,i?)«-(«)-'??([AO,[s,e]), 

that proves if^ . I 



Proposition 3 There is an identification of the orthogonal complement of the vector x Cz Ti. with the 
tangent space to the U{TC)-orbit through ^ = fJ-{x) in u*{Ti.). For y,y' G Ti. orthogonal to x with respect 
to the Hermitian product, the vectors {lJ-*)x{y), {p-*)x{'y') o-f£ tangent to the orbit through ^ and 

'^li{i^*Uy),ii^*Uy')) = 9iy,y'), (23) 
vli,ii^*Uy')) = ^iy,y'), (24) 
"-MxXi^^Uy)) = ifi*)AJy)- (25) 



9 



Proof. Since 

ifi.Uy) = Py =\y){^\ + \^){yl 

II 11^ 

can be written as Py — [Ay, where Ay is a Hermitian operator such that Ax = iy and Ay = —ij^x, 
the operators Py,Py,, viewed as vectors in u*{H), are tangent to the orbit through ^. Then, due to 

<-liPy^Py') - ^Tr(P;oP;,) = ^Tr(I|.f.|y)(j/'|+(y,2/%.|x)(.|) 

^ {\\xr{{y',y)n + {y,y')n))^my,y')n)^g{y,y')- 



2 X 



To prove we use (tTzjl : 



^l(Py^ Py') = -(^' lA' ^.'])n*(«) = - ^Tr o [A„ Ay, 



2 



(a;, [Aj^, = ((^ys;, Ay>x)-H - {Ay,x, Ayx)-H) 



FinaUy, foUows directly from uj{y', Jy) = g{y',y) and I 

The above theorem says that the Kahler manifold {Vl {7i),^J^, cr'', tj^) comes from a sort of a 'Kahler 
reduction' of the original linear Kahler manifold (TYr, J,g,Lu). In particular, the symplectic manifold 
V^.iH) is the symplectic reduction of (TYrjU;) with respect to the isotropic submanifold Sr = {x ^ H : 
\\x\\ = r}. The characteristic foliation of consists of orbits of the group = {z ^ C : \z\ = 1} 
acting on Ti by multiplication. The fundamental vector field of this action is 

-U = pkdq, - qkdp, 

which is simultaneously a Killing vector field for the Riemannian metric g. Therefore g induces a 
Riemannian metric on I?^(7i), etc. 



3 Smooth manifold structure on V^{7i) 

Recall that the space of non- negatively defined operators from gl{TC), i.e. of those p G gli^H) which 
can be written in the form p = T^T for a certain T e glCH), we denote by V{'H). It is a cone as 
being invariant with respect to the homoteties by A with A > 0. The set of density states ViTi) 
is distinguished in the cone V{Ti.) by the equation Tr(p) = 1, so we will regard V{H) and 'D{T-L) as 
embedded in u*{H). 

The space V{Ti.) is a convex set in the afhne hyperplane in u*{Ti), determined by the equation 
Tr(T) = 1. The tangent spaces to this afHne hyperplane are therefore canonically identified with the 
space of Hermitian operators with trace 0. It is known that the set of extreme points of 'D(Ti.) coincides 
with the set 2?^(7i) of pure states, i.e. the set of one-dimensional orthogonal projectors \ x){x \ (see 
Corollary O . Hence every element of 2?(7i) is a convex combination of points from 'D^{H). The 
space T>^{TL) of pure states can be identified with the complex projective space PTL ~ CP"~^ via the 
projection H \ {0} 9 x i— > | x) (ccj S T>^{H) which identifies the points of the orbits of the C \ {0}-group 
action by complex homoteties. We have already seen that 'D^{H) is canonically a Kahler manifold. 
This will be the starting point for the study of geometry of the set 'D{H) of all density states. 

The (co)adjoint action of the group U{H) in u*{TC) induces its action on the positive cone ViTi.) and 
on the space of density states. This action is transitive on pure states but it is no longer transitive on 
subsets V''(n),k>l, where V''(n) = V{n)f\V^{n) and V^{n) consists of non-negative operators of 
rank k. The rank is understood clearly as the rank of the corresponding operator (or matrix, if a basis 
in TL is chosen). The intersection of 'D{'H) with any Weyl chamber in a Cartan subalgebra in u*{7i) is 
an {n — l)-dimensional simplex, while the intersection of ^^{TL) is the (fc — l)-skeleton of this simplex. 
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However, the dimension of the orbit may vary even for points from a chosen V^{rL) if fc > 1. Thus, the 
set of density states is a union of smooth manifolds ~ orbits of U{T-L) - but the differentiable structure 
of the stratum D^iTC) is a priori not clear (for fc > 1), since the decomposition into orbits is not a 
regular foliation, i.e. V^iTL) is the union of a family of various submanifolds of u*{Ti.) which differ even 
by dimensions. By the differential structure we mean here the differential structure inherited from 
u*{Ti.), so that the smooth curves in 'D(Ti.) and hence the tangent spaces are uniquely defined. 

Our aim in this section is to understand this differential structure. Of course, the interior of 
namely X'"(7i), is an open subset, so a submanifold, in the affine subspace of trace=l Hermitian 
operators and the real question is only the boundary, consisting of those density states p for which 
det(p) = 0. The best situation would be if the boundary were a submanifold, but this is not true in 
dimensions n > 2 as we will show later. The stratification into [/(7i)-orbits is too small, since, as it 
will appear later, the subsets V^{TL) are coarser submanifolds in u*{TC). We will show also that the 
stratification by rank is the maximal one in the sense that the vectors tangent to 2?(7i) at p G I?'^(7Y) 
must be tangent to 'D^{'H) itself, so the largest it* (7i)-submanifold through p G 'D^{'H) contained in 
V{H) is V^{H). 

We start with fixing an orthonormal basis in H. which allows us to identify u*{TC) with the space 
u* (n) of Hermitian n x n-matrices which is canonically an n^-dimensional real manifold with respect 
to the identification 

u*in) 3 (a,,) ^ ((a,,)?, (a.,),<,) G R" x C«("-i)/^ 

By V{n) we denote the space of non-negatively defined matrices from u*{n), by V^{n) the subset of 
rank fc matrices from V{n), etc. Let us denote by 'Pj{n) the set of matrices A — {aij)^,j^i G Vin) being 
of rank fc and such that the minor det[(ars)r,sej] associated with a set of indices J = {ji, . . . ,ik} C 
{1, . . . , 7i} is non- vanishing.^ The next lemma shows that any matrix from Vj{n) can be reconstructed 
from its rows (or columns, since it is Hermitian) indexed by J. 

Lemma 2 Let A — {aij)'2j^i G Vj{n), so that the matrix {ars)r.seJ has the inverse {a'^^)r,sej ■ Then 
the matrix A is uniquely determined by {{aij) : i E J, j = 1, . . . ,n} according to the formula 

ai-j = ^ airU^'^aj;. (26) 

Proof.- The matrix A being non-negatively defined is of the form T^T for certain n x n-matrix T, so 
that Oij is the Hermitian product {ai,aj) of columns of T with respect to the standard Hermitian 
product (O. The matrix T is not uniquely determined. However, the fact that A is of rank fc with 
the non- vanishing minor associated with J means that the columns aj , j G J are linearly independent 
and span the rest of the columns of T. But the Hermitian product on the subspace in C" spanned by 
{ttj : j G J} is given by the formula 

{x,y)c^- ^ ^ (x,ar)c"a'''*(as,y)c", (27) 

where (a'^'')r^seJ is the inverse of the matrix ((a^, as)c" )»-,seJ- The proof of (|27|l is immediate, since 
the r.h.s. of H27() is C-linear with respect to y, anti-linear with respect to x and equals (ai,aj)c" for 
X — ai, y — Uj, i,j G J, by definition. Since aij = {ai,aj)c", we get the formula (|26|l directly from 
(EZ|l- ■ 

Remark. It is worth noticing that the formula (|27|l is similar to the one describing the Dirac bracket 
on constraint manifolds induced by second class constraints. 

For J as above define a linear map 

$j : u*{n) u*(fc) X C^"-*^)*^ ~ R*^ x c'^^^^k-k'-k)/2 ^ ^2nk-k^ 

^The set J is not to be confused in the following with the complex structure denoted accidentally by the same letter. 
From the context, however, the notion of J is always obvious. 
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by 

*j((a»i)"i=i) = ((a*i)ij6J, {ars)r<^j,sej))- (28) 
In particular, if we work with the principal minor, i.e. J = {!,..., A;}, then $j associates with a 
Hermitian matrix its first k columns with removed, say, upper-triangular part which is irrelevant due 
to hermicity or, equivalently, its first k rows with removed lower-triangular part. 
For A E u*{n) by ^j^a we denote the map ^j^a{X) = - 

$j,A((a;jj)"j=i) = {{xtj - a^j).i^j(zj, {xrs - ars)r(^j,sej))- (29) 

With some abuse of notation, its restriction to V^{n) we will denote by the same symbol. It is clear 
from the above Lemma that the map $j is continuous and injective on Vj{n). Thus, for A G 'Pj{n), 
also the map is continuous and injective on 'Pj(n). 
Conversely, every point 

{{yij)i,]eJ, {yrs)r(^.J,sej)) 

of u*{k) X C^""'')'' ~ R^"'""'"^, sufficiently close to 0, is the value $j_^(X) for a certain X e Vjin). 
Indeed, adding a small Hermitian matrix to (aij)ij^j will not change its invertibility. Hence we have 
to reconstruct X out of $j(X), i.e out of the columns (and rows, since X should be Hermitian) with 
indices belonging to J and knowing that {xij)ij^j has an inverse, say, {x^^)r,sej- Here Xij = aij +yij 
for j e J. An obvious choice is the formula i.e. 

Xij — ^ ^ Xij'X Xjs- 

The only thing to be checked is that X ~ {xij)fj^i defined in this way is non-negatively defined 
and of rank k. Assume, for simplicity of notation, that J = {1, . . . , k}. First, we can find vectors 
Pi, . . . ,l3k & such that 

Xij = {Pi,P])c'' (30) 

for i, j — 1, . . . , k. This can be done up to a unitary transformation. For example, (3i can be columns 

of the matrix ^Jixij^t^i- Then, we find (this time unique) vectors Pk+i, ■ ■ ■ , Pn G satisfying the 

conditions Xij = {Pi, Pj)c'' i * = k+1, . . . ,n, j = 1, . . . ,k. It is easy to see now that, due to the formula 
(|27|l . we have for all i,j = 1, . . . ,n. This immediately implies that X is non-negatively defined 
and of rank k. Moreover, since 

Xij = ^ {a^r +yir)al^{yj; +a~), (31) 

where {ay'^)r,seJ is the inverse of the matrix (a^s +yrs)r,sej, the matrix elements Xij rationally depend 
on ymi, so that is smooth, thus also regular, as a function from a neighbourhhod of in R,2nfc-fc^ 
into u*(n), so V^{n) is a submanifold in u*{n). To see the image of the differential of ^■'^ 
tangent space TAV^in), let us consider the linear (with respect to y) part (wy ) of the r.h.s. of (|31|l : 

E/ rs rm Is . rs \ /oo\ 

(j/irO Qjs - QirU y^ia ttjs + ttirtt y^s)- (o2) 

r,s£j 

To see this better, let us change the orthogonal basis of C" for such that J = {1, . . . , A:} and A is 
diagonal, an = A^, Ai = for i > k. Then one can easily find that l|32|l takes the form 

(33) 

This means that in the image are arbitrary Hermitian matrices V = (I'yOfj^i such that Vij — for 
i,j > k, that can be written in a coordinate-free way as {Vx,y)c'^ = for all x,y £ Ker(^). Note 
that the manifold V''{H) is connected. Indeed, it consists of connected orbits of the group U{H) 
which meet a Weyl chamber as the (fc — l)-dimensional skeleton of a simplex. However, the connected 
components of this skeleton are identified by the action of the Weyl group, so they form topologically 
a (fc — l)-dimensional simplex which is obviously connected. Therefore we have proved the following. 
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Theorem 1 Let A e Vjin). Then the map $xa : V^ln) R2"fc-fc defined by H^j is a local 
homeomorphism from a neighbourhood of A in V^in) onto a neighbourhood of in u*{k) x C^""*^^*^ ~ 
p^2nfc-/c ^ Moreover, the collection of the maps 

^TlA ■ ^J,A V^{n) C u*{n) defined on sufficiently 
small neighbourhoods Wj.a of by the formula \S1]) constitutes a smooth manifold structure on V^{n) 
which makes it into a smooth and connected submanifold of u*{n). The tangent space Tyi'P'^(n) , viewed 
as a subspace of u*{n) consists of matrices V G u*{n) satisfying (Vx,y)c" ~ for all x,y £ Ker{A). 

Remark. In section 5 we obtain the manifold structure on V^{n) much simpler as the structure of 
an G'i(n, C)-orbit. But we find that Lemma El and Theorem Q are of some interest per se providing 
explicit coordinate systems. 

The next theorem shows that smooth curves in u*{n) which lay in T'{n) cannot cross V^ln) transver- 
sally, i.e. V^{n) is in a sense an edge for V^^^{n) if fc < n — 1. 

Theorem 2 Let 7 : R — > u*(n) be a smooth curve in the space of Hermitian matrices which lies 
entirely in V{n). Then 7 is tangent to the stratum V^{n) it belongs, i.e. j{t) G V''{n) implies 
7WeT^(t)P'=(n). 

Proof - Of course, it is enough to prove the above for an arbitrary t gH, say, t = 0. Assume therefore 
that A = 7(0) e r''{n). Take x £ Ker(A). Since 

for At > 0, we have {■^{0)x, x) > 0. Taking in turn At < we see in a similar way that (7(0)x, x) < 0, 
so 

(7(0)x,x)=0. (34) 

By polarization of (|34|l we get 

(7(0)a;,y) + (7(0)y,x) =0 (35) 
for all x,y Ker(v4). But 7(0) is Hermitian, so 

{j{0)y,x) = (y, 7(0)2:) 

and H35() yields that the real part ^{{j(Q)x, y)) is for all x,y £ Ker(A). On the other hand, the kernel 
of >1 is a complex subspace and 

5R((7(0)x,i.y)) =3((7(0)x,y)) 

so 

(7(0)a:,y)=0 (36) 
for all x,y G Ker(A). But, according to Theorem^] H36() means that 7(0) £ T^7''^(n). I 



4 Smooth stratification of density states 

The set 2?(7i) of density states on H is the intersection of the cone V{H) with the affine subspace 
{A £ u*{H) : Tt{A) = 1} or, in other words, it is the level set of the function Tr : ViH) -> R 
corresponding to the value 1. Since Tr(tp) = iTr(p) and V'' is invariant with respect to homoteties 
with positive t, it is clear that Tr is a regular function on each V''{'H), so that V^iTi) is canonically a 
smooth manifold. Since topologically V'^{n) ~ V^{n) x R, the manifolds V'^{n) are connected. All 
these observations together with Theorems ^ a-nd |21 can be summarized in the following. 

Theorem 3 The spaces T>^{Ti.) of density states of rank k, k — 1, . . . ,n, are smooth and connected 
submanifolds in u*{H) of (real) dimension 2nk — k^ — 1. The tangent space TpV'^{'H) is characterized 
as the space of those Hermitian operators T of trace which satisfy {Tx,y) = for all x,y € Ker(p). 
Moreover, the stratification into submanifolds T>^{T-C) is maximal in the sense that every smooth curve 
in u*{H), which lies entirely in T>{Ti,), at every point is tangent to the strata T>^{'H) to which it actually 
belongs. 
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Corollary 2 The boundary O'DiTi) — Ufc<n ■^'^(^) '^f density states is not a smooth sub- 

manifold ofu*(Ti.) if n = dimH > 2. 

Proof - If n > 2 then the boundary d'D{H) has at least two different strata and the vectors orthogonal 
to, say, the stratum 'D^{H) of pure states are not tangent to dT>{Ti.). But the dimension of V^lTL) is 
smaller than the topological dimension of d'D{H). I 

Remark. It is well known that for n = 2 the convex set of density states is affinely equivalent to the 
three-dimensional ball and its boundary - to the two-dimensional sphere, so it is a smooth manifold. 

The last problem concerning the geometry of density states we will consider is the question of afhne 
parts of the manifolds T>''{H). It is motivated by the fact that the set T>^{H) of pure states is exactly 
the set of extremal elements of X'(7i), so it does not contain intervals, but the other strata T>''{H) with 
k > 1 must do as shows the following theorem. Recall that a non-empty closed convex subset Kq of 
a closed convex set K is called a face (or extremal subset) of K if any closed segment in K with an 
interior point in Kq lies entirely in Kq; a point x is called an extreme point of K if the set {x} is a 
face of K. 

Theorem 4 If p e V^{H) then the affine space in u*{TL) which is tangent to Ty^iJ-L) at p intersects 
'D{'H) along a (fc^ — 1)- dimensional convex body which is affinely equivalent to the set of density 

states in dimension k. This convex body is exactly the face ofT>{Ti.) at p. In other words, the face of 
T>{TC) at p ^ T>^{l-C) is affinely equivalent to T>{k). 

Proof.- Let us take coordinates in u*{Ti,), i.e. let us chose an orthonormal basis in H, in which p is 
represented by a diagonal matrix {pij), pij — S^Xt, where Xt — for i > k. According to the form of 
TpV^{TL), matrices {xij) which belong io p + Tp'D^{T-L) have entries Xij with i,j > k equal to 0. If 
they belong as well to T>{Ti.), also Xij — if i > k or j > k. Indeed, since Xij = {ai,aj) for certain 
vectors Zi 6 C", we have xu = = 0, so at = 0, for i > k, and further Xij = {at, aj) = if i > fc 

or j > k. In other words, the only non-zero part of X is the block (a;ij)f ^-^i which is therefore an 
element of T){k). Conversely, every matrix X with such a block form belongs simultaneously to V{H) 
and, since {X - p)„ = for i,j > k,to p-\- TpV^iU). To see that {p + TpV^{'H)) n V{rL) is exactly 
the face of T>(T-C) at p, consider a segment in T>{Tl) for which p is an interior point. The open segment 
is clearly a smooth curve in T>{7i), so, in view of Theorem|ni it is tangent to T>^{TL) at p, thus belongs 
entirely top -f-TpI3'=(7^). I 

Corollary 3 Extremal points of T>{TC) are exactly pure states. 



Let us mention that a major part of what has been said about the differential structure of the space 
V^{Ti.) of rank-fc positive operators can be repeated for the space of all rank-fc Hermitian operators. 
Denote by u^^ f. (Ti.) the set of those Hermitian operators ^ whose spectrum contains k^ positive and 
fc_ negative eigenvalues (counted with multiplicities) , respectively. Thus the rank of ^ is fc = fc+ -I- fc_ 
and pfc(n)=< o(n). 

Fixing an orthogonal basis in TC will identify u^^ ^.^ (Ti) with the space u^^ (n) oi n x n Her- 
mitian matrices of rank fc with the corresponding spectrum. Denote by I?^^ the diagonal matrix 
<im(7(l, . . . , 1, — 1, . . . , — 1, 0, . . . , 0) with 1 coming fc+-times and —1 coming fc_-times. Denote by 
(•, •)fc+,fc„ the 'semiHermitian' product in C" represented by -D^^: 



5 Geometry of u*{l-L) 



k++k- 




(37) 



It is easy to see the following. 



14 



Proposition 4 Any ^ = {(^ij) G k "^"^^ written in the form ^ = T'^d'^^T for eertain T G 
GL{n, C). /n other words the entries of the matrix ^ are semiHermitian products Oij ~ {ai, Oij)k^,k- , 
where ai denotes the ith column of T . 

Proof. We can diagonalize ^ by means of an unitary matrix U , 

U^U^ = diag{Xi, . . . , A„), 

wliere Ai > • • • > A„, so Ai, . . . , A^^ > and Afe++i, . . . , Xk^+k_ < 0. Hence ^ = T^l+T for T ^ CU 
witli 



C = diag {y\Xi\, ^\Xk++k^\, 1 ■ • • , l) 



Now, we can reformulate Lemma |21 for u^^ fe_ (^i) instead of V^{n). Tlie proof is essentially the 
same with the difference that we use the semiHermitian product (•, ■)fc+,fc_ in C" instead of (•, •)c"- 

Lemma 3 Let^^ = {aij)fj^-^ £ u^^ {n). Assume that the matrix {ars)r,seJ has the inverse {a^^)r,seJ 
for certain k — (fc+ + k^)-element subset J = {ji, . . . , jfc} C {1, . . . , n}. Then the matrix ^ is uniquely 
determined by {(ay) : i G J, j = 1, . . . , n} according to the formula 

ai-j = ^ aira^'^aJJ. (38) 

One can now prove that (7i) are submanifolds of u*(J-C) in completely parallel way to the case 

of V^iTi). However, Proposition^] suggest an easier (although less constructive) way to do it. Namely, 
we can see ^_ {7i) as an orbit of a natural GL{'H) action on u*{TL). 

Theorem 5 The family 

{wfc+,fc_ (^) : k+, k^>0,k = k+ + k^ <n} (39) 
of subsets of u*{Ti.) is exactly the family of orbits of the smooth action of the group GL{Ti.) given by 

GL{n) X u*{n) 3 (r,o ^ t^t^ e u*{h). (40) 

In particular, every u'^^ (Ti) is a connected submanifold ofu*{Ti) and the tangent space to u^.^ (Ti) 
at ^ is characterized by 

B e T^ul^^k_{H)^yx,ye KeriO [{Bx,y)n = 0]. (41) 
Moreover, the following are equivalent: 

(1) u*f.^ fc_ ('^) intersects P{H); 

(2) u^^ fc_ ('^) is contained in V{Ti,); 

(3) fc_ = 0; 

(4) ul^ ,_{n)=V'^{n),k = k+ + k,. 

Proof. The proof that (|40|l is a group smooth action is straightforward. Proposition 0] shows that 
ul (H) is contained in the GL(7i)-orbit of d'^^. 

On the other hand, although the spectrum is not fixed on every G'i(7i)-orbit, the number of positive 
and the number of negative eigenvalues (counted with multiplicities) are fixed along the orbit. Indeed, 
if {x,^x)-H > (resp. {x,S,x)-h < 0) for a; in a fc+-dimensional (resp. fc_-dimensional) linear subspace 
V+ (resp. VI), then {x,TST'^x)n = {T^x,^T'<x)n > (resp. {x,T^T^x)h = {T'<x,^T'<x)n < 0) for x 
in a fc_|_-dimensional (resp. fc_ -dimensional) linear subspace {T'')^^{V+) (resp. {T'')~^{V-)). 
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The corresponding infinitesimal action of w £ gl{'H) is f i-^ + and the operators — w^ + ^w^ 
clearly satisfy {x,^yy)T-c = for all x,y G Ker(^). Conversely, if for certain B £ u*{Ti) we have 
{Bx,y)-H = for all x,y <E Ker(f), then B can be written in the form v£, + ^w^. To see this, consider 
the splitting H = Pi V2, where V2 = Ker(^) and Vi = . According to this splitting ^ can be 
written in the operator matrix form 

■ a 





where ^1 is Hermitian and invertible. Similarly, B has the form 

B -- 



Bii B12 
B21 

where bI^ — Bn and, in the obvious sense, B21 = B\2. Now, it is easy to see that _B = + £^v\ 
where 

i?2ier' 

that proves H41|) . 

Finally, if {TL) intersects V{Ti.), then it contains an element with non-negative spectrum. But 

the signs of the elements of the spectrum are constant along a Gi(7i)-orbit which means that k- = 

and ui^^,_{n)^rHn)crin). i 

Note that the fundamental vector fields a(^) — —a£,—£,a^ of the G-L(?i)-action satisfy the commutation 
rules \a,b]yf — [a, 6]_ 

The next results shows that the foliation into submanifolds uj!^ {H) can be obtained directly from 
tensors A and R. We know already that the (generalized) distribution Da induced by A is generated by 
vector fields Aa{£,) — #A^(A) = [A,S,] and the (generahzed distribution Dji induced by R is generated 
by vector fields Ra{0 = = + - The following is straightforward. 

Theorem 6 The family {Aa, Ra ■ A G u*{'H)} of linear vector fields on u*{Ti.) is the family of 
fundamental vector fields of the GL{7i)- action: 

AAiO = = (42) 

Ra{0 = A^ + ^A^A^ + ^A^^-A{0- (43) 

In particular, 

[Aa, AB]vf — A[A,B], [Ra, RB]vf — A[A,B], [Ra, AB]vf — R[A,B], (44) 

so the (generalized) distribution induced by jointly by the tensors A and R is completely integrable and 
fc (^) '^''6 maximal integrate submanifolds. 

Corollary 4 The generalized distributions Dgi — + D\, D\ and Dq = Dfjf]D\ on u*[TC) are 
involutive and can be integrated to generalized foliations Tgi, Tt^, and Tq, respectively. The leaves of 
the foliation Tgi are the orbits of the GL{Ti.) action ^ f— > T^T\ the leaves of J- a are the orbits of the 
U{TL)- action. 

Denote by J and the (1, l)-tensors on u*{'H), viewed as a vector bundle morphism induced by 
the contravariant tensors A and R, respectively, 

J,R : Tu*{n) ~*Tu*{n), 
J^A) = [A^] = A^{A), 
n^A) = [A,^U^R^{A), 

where A e u*{H) ~ T^u*{H). The image of J is Da and the image of TZ is Dr. 
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Lemma 4 The tensors J and TZ commute and 

o {A) = 7^^ o (A) = [A, e] ■ (45) 
Proof. We have ^ 

But, as easily seen, ^ ^ 

+ = = [[A,C],C]+ - o (46) 

I 

Recall that J7(7i)-orbits O, i.e. the orbits with respect to the action of the subgroup U{Ti,) C 
GL{Ti)^ carry canonical symplectic structures rf^ . The symplectic structures rP is [/(7i)-invariant, 
i.e. {O^rj'^) is a homogeneous symplectic manifold. We will show that this symplectic structure is a 
part of a canonical Kahler structure. We know already this structure for the orbits V^{H). 

Recall also that on u*{Ti.) we have the Riemannian metric induced by the scalar product {A, = 
iTr(AB) on u*{H). 

Theorem 7 (a) The image of is T^O and Ker{J^) is the orthogonal complement of T^O. 

(b) is a self-adjoint (with respect to {■,-)u* ) md negatively defined operator on T^O. 

(c) The (1, l)-tensor J on u*{Ti.) defined by 

JdA) = {-{JOhoY' ^dA) (47) 
induces an U{7i) -invariant complex structure J on every orbit O. 

(d) The tensor 

^f{AB)^r^f{A,J^{B)) (48) 
is an U {Ti)- invariant Riemannian metric on O and 

^fiJ^iA),B)^fjfiA,B). (49) 

In particular, (O, ry'^, 7'^) is a homogeneous Kahler manifold. Moreover, if £, G u*{7i) is a 
projector and ^ G O, then — and 7'-' (A, B) — (A, B)^* . 

Remark. The tensor J is canonically and globally defined. It is however not smooth as a tensor 
field on u*(fK). It is smooth on the open-dense subset of regular elements and, of course, on every 
C/(7Y)-orbit separately. 

Proof. 

(a) The vector fields A^(^) = \A,£\ = J^£^{A) are fundamental vector fields of the [/(7Y)-action, so 

T^O is the image of J'^. Moreover, the invariance of the Riemannian metric {A, B)u* , 

{J^{A),B)^, = {[A,^],B)u. = -{AJ^{B))u', (50) 

implies that 

SeKcr(Jj) ^ B^J^{u*{H)). 

(b) The identity H5U|) means that = — J7^, where is the adjoint operator to with respect to 

the scalar product (A, . Consequently, 

{Jir^Jl (51) 

Moreover, is negatively defined on T^O, since 

(j|(A),A)„. = ([[A,e],e],^)«- =-([A^],[A^])«* <o, 

for eT^O, [AC] ^0. 
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(c) The tensor J is clearly J7(7^)-invariant: 

Jv(,uAUAU^) - \UA\j\VilJ^\ = U{A,i\\J^ = lJ{Je^{A))U\ (52) 

so J7(7i)-invariant is the tensor ^— J'^^ ^ and its composition J . The tensor J' defines an 
almost complex structure on every orbit O, since 



To show that this almost complex structure is integrable, it is sufiicient to show that the dis- 
tribution M in the complexified tangent bundle TO ® C which consists of i-eigenvectors of 
(complexified) J is involutive. Since J ^ and therefore A/", is invariant, it is sufficient to check it 
at one point, say <^ S O with respect to the complexified Lie algebra gliH) = u* (Ti) C equipped 
with the bracket [a, b] = l[ab — ba]. 



Let — Kj, . . . , — where ki, . . . , Km > 0, be the eigenvalues of {J^)\T(:0 counted with multi- 
plicities. The complexified J^, which with some abuse of notation we will denote by the same 
symbol, has therefore eigenvalues iiKk with eigenvectors a^, k = 1, . . . , m and ^/^(a^) = i^^fc • 
Thus A/j is spanned by the vectors a'^ ,k = 1, . . . , m, i.e. eigenvectors of J^, J^{al) — inka'^ 
with positive Kk- This space is clearly a Lie subalgebra in glili.), since 

the vector [a^, a^], if non-zero, is again an eigenvector of corresponding to a 'positive' eigen- 
value i{Kk + Kl). 

(d) The tensor 

jfiA,B)^r^f{A,J^m 

is clearly [/(7i)-invariant. From H5()(l and H51(l it follows that J'^ — —J^- Since J and J clearly 
commute, J^d^'^l) = in view of JHIl, 

^'^{{AA.Jmm = ([Ae],J?(S))u*(«)-(-Jc([A^]),S)u-(w) (53) 
= -r7f(Jc([A,e]),[S,e]). 

This immediately implies that is symmetric and proves H49|l . But (|17|l implies also that 

7f([Ae],[Ae]) = ryf([Ae],Jc([Ae]))-([Ae],Jc(^))«n«) (54) 
= (A-JjJ?(v4))„.(w). 



But 



is a positive operator, so 



7^^([A^],[A^])>o 



for [AC] ^0. 

Finally, if ^ is a projector, ^ — ^, then (cf. (|^ ') 

Jj2([AC]) = -[AC], 

so = and (cf. Q) 

7f ([A,C],[i?,^]) - ([AC], J«(i?))„.(«) - {{A,i\ABX\)u'iu). 
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We have some similar results for the tensor TZ which however are not completely analogous, since 
the distribution Dn is not globally integrable. The proofs are analogous, so we omit them. 

Theorem 8 (a) The image -Dfl(C) of TZ^ is the orthogonal complement of Ker{TZ^). 

(h) TZ^ is a self-adjoint (with respect to (■, ) and positively defined operator on 

(c) The (1, 1) -tensor TZ on u*(Ji.) defined by 

7^e(A) = |(7^e)p„(c)r'°^^(^) (55) 

satisfies TZ^ = TZ. 

Corollary 5 The distribution Dq is the image of J^^oTZ^ = TZ^oJ'^. In other words, i?o(0 = {[^i^^] • 
A G u*{Ti)}. Moreover, the foliation J-q is U{Ti)- invariant, J -invariant and TZ-invariant, so that J 
andTZ induce on leaves ofJ-Q a complex and a product structure, respectively. The leaves of the foliation 
J-Q are also canonically symplectic manifolds with symplectic structures being restrictions of symplectic 
structures on the leaves of J-\, so the leaves of J-q are Kahler submanifolds of the U{Ti)-orbits in 
u*{n). 

Proof. The image of o TZ^ — TZ^ o is clearly contained in Dq. Conversely, let B G Dq{^) = 
^Apl^R.- According to H46(l . -Do(C) is invariant with respect to both: and 7?.^ and and 7^^ 
are injective, thus surjective, on Dq{^). The distribution Dq is therefore generated by vector fields 
Xa{£,) = [A,^'^]. It is a matter of simple calculations to show that these vector fields commute with 
the fundamental vector fields of the U{Ti) as [XAT^B\vf — A[B,yi] that shows U{TC) invariance 
of Dq. One can also easily seen that the restrictions of ^'■^ to the leaves of JFq contained in O are 
non-degenerate. It follows also directly from the explicit calculations we present below. I 

Let us explain the above theorem in local coordinates, i.e. for the case of matrices. Suppose that 
^ = diag{Xi, . . . A„) € u*(n) is a diagonal matrix. For simplicity, it is better to start already with the 
complexified structures, i.e. with gl{n) = u*{n) (g) C equipped with the bracket [a, b] = j{ab — ba) and 
the Hermitian product {a,b)gi = iTr(a^6), so that u*[n) is a real Lie subalgebra in gl{n) with the 
induced scalar product. Let Ef be the matrix whose the only non-zero entry is 1 at fcth row and Zth 
column. We have 



and 

so that 
and 

In particular, 
Consequently, 
and 
so that 



{Ef,E:),i^\{5l5[), 


(56) 


[E^,E:]^~t{SlE^~6':El), 


(57) 


[EtE:U^6lE'; + S',El, 


(58) 


J^{Ef) = [El£]=^{Xk-\i)Ef. 


(59) 


TZ^{E^) = [Eli]+^{\k + \i)El 


(60) 


J^oU^iEh = [Ete] = ^i>^l - >^f)Et 


(61) 


Jl{E^) = -{X^^\ifEt 


(62) 


nl{Ef)^{\k + XifEt 


(63) 


J^{Ef)=i-sgn{Xk-\i)El 


(64) 
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and 

n^{E^) = sgn{Xk + Xi)E^. (65) 

The (complexified) tangent space T^O ® C is spanned by those Ej' for which — A; 7^ 0, the space 
-Dfl(C) <8) C is spanned by those E^ for which Afc + A; 7^ 0, the space -Do(C) C is spanned by those 
E^ for which Xl — Xf ^ 0, and the distribution TV mentioned in the proof of the theorem is spanned 
by El for which A/j — A; > 0. The complexified symplectic form reads 

vf{i{Xk - Xi)Eli{Xr - X,)El) = (i(Afe - Xi)ElEl),i = -i{Xu - Xi)\{5l5^,), 



I.e. 



and the complexified Riemannian form 

^f{El El) = rjfiEt JdK)) = 2|a/-A,| ^'^^'^^^- ^^^^ 
As a basis in u*{n) let us take 

A'l = Ef + El k<l, Bf=iEf-iEl k<l. (68) 
It is easy to see that this is an orthonormal basis and that 

J^{Af) = sgn{Xk - Xi)Bl J^{B^) = sgn{Xi - Xk)Af. (69) 

and 

TZ^iAf) = sgniXk + Xi)Af, J^{Bf) = sgn{Xi + Xk)Bf. (70) 

Moreover 



nf{Bl Al) = j^^y nf{Bl Bl) = vfiAf, a:) = 0, A, - A,, A. - A, ^ 

and 



and 



(71) 



7f(Bf,AD=0, 7f(Bf,B;)=7f(A',^D = ^^^> A,-A,,A.-A,^0. (72) 



In other words 

^f= E I, \ , (dbf ^d6f + daf ^daf), (74) 
where 

are coordinates on u* (n) such that Bf^ = df,k , Af = d^k . The reduction of the symplectic form rf^ to 
the leaves of the foliation 

is clearly non-degenerate and constitutes, together with the reduced Riemannian structure 

(7f)|^o= E n-^(d6f®d6f + daf®daf), (76) 
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a Kahler structure. 



Remark. Of course, when ^ is a projector, then Afc — 1,0, so Afc — A; ^ =^ |Afc — Aj | = 1 and 
reduces to the canonical scalar product. Note also that the leaves of T\ and JFg through ^ coincide, 
except for the rare case when there are A, A' 7^ in the spectrum of ^ such that A+A' — 0. In particular, 
the foliations J-\ and J-q coincide when reduced to the subset ViTL) of non- negative operators or to 
the set 'D[T-L) of density states. On such leaves the product structure TZ is trivially the identity. 

6 Composite systems and separability 

Suppose now that our Hilbert space has a fixed decomposition into the tensor product of two Hilbert 
spaces Ti = Ti.^ (S> Ti.^- This additional input is crucial in studying composite quantum systems and it 
has a great impact on the geometrical structures we have considered. The rest of this paper will be 
devoted to related problems. 

Observe first that the tensor product map 

<^ : X ^ H = (77) 

associates the product of rays with a ray, so it induces a canonical imbedding on the level of complex 
projective spaces 

Seg : PH} X PH^ PH ^ P{n^ ® U^), (78) 

{\x^){xWx'){x^\) ^ (E)x^){x^ (g)x^\ . (79) 

This imbedding of product of complex projective spaces into the projective space of the tensor product 
is called in the literature the Segre imbedding The elements of the image Seg{PH^ x PTi^) in 

PTi = P{'H} (g) Ti?) are called separable pure states (with respect to the decomposition Ti. — Ti} ® Ti^). 

The Segre imbedding is related to the (external) tensor product of the basic representations of 
the unitary groups U{Ti}) and U{Ti?), i.e. with the representation of the direct product group in 

{p^ ® p^){x^ ® x^) = p^{x^)® p^{x^). 

Note that p^ ® p^ is unitary, since the Hermitian product in TL is related to the Hermitian products in 
V} and by 

{x^ <»x'^,y^ (»y'^)n = {x\y^)w ■ {x\y^)n^. (80) 
The above group imbedding which, with some abuse of notation, we will denote by 

Seg: U{H^) x U{n^) U{n), 

gives rise to the corresponding imbedding of Lie algebras 

Seg : u{H^) x u{n^) u{n), 

or, by our identification, of their duals 

Seg: u*{H^)xu*{H^)^u*{H). (81) 

The original Segre imbedding is just the latter map reduced to pure states. In fact, a more general 
result holds true. 

Proposition 5 The imbedding fUJ) maps V'^iH^) x P'{n^) into V^^rO- ®n'^) and V^{V}) x V^Ti"^) 
into V^\n^®T-e). 
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Proof. Let us take e V^{U^) and e V^V?). Using bases of eigenvectors {e\) of and (e|) of 

to construct a basis (e\ (gj e'j) of eigenvectors of (8) one easily sees that the elements of the 
spectrum of A^ (g) A"^ (counted with multiphcities) are A,A^-, where A^{el) = Xic] and ^^(e|) = Aje|, 
so that A^ iSi A^ = Seg{A^ , A^) is non-negatively defined and has rank kl. If A^, have trace 1, i.e. 

A, = 1 and E,- a;- = 1, then E^,,- A^A^- = A^ • E,- = 1- ■ 

Let us denote the image Seg{V''{n^) x V^iU"^)) by S^'\n^ ®n^), the set ® H^) of separable 

pure states simply by S^{H^ (8> W^), and the convex hull 

conv(Seg {V{n^) x ^(W^))) 

of the subset Seg (P(H^) x V{H'^)) of all separable states in by 7^^). The states from 

i.e. those which arc not separable, arc called entangled states. 

Proposition 6 The convex set S{Ti} ®H^) of separable states is the convex hull of the set S^ili} ^H^) 
of separable pure states and iji} ^TH?) is exactly the set of extremal points ofS{H^®'H'^). Moreover, 
S^CH^ (8) W^), thus S{H^ <S) 'h?), is invariant with respect to the canonical Uiji}) x IKJi})- action on 

(Ti, T2)A = (Ti O T2) o A o (Ti ® T2)t. 

Proof Let us start with showing that the convex hull of S^iW^ O 'H?) contains Seg {V{'hO-) x Vi'hP)) 
thus equals S{'H^ ® Ti?). Indeed T>^(W) is the set of extreme points of 2?(W), i = 1,2, so that any 
A* € X>(W*) is a convex combination A' = tlpl of elements pi G V^{W), i = 1,2. Hence, A^ (g) A^ is 
the convex combination 

A^ ®A^ =Y,t\tl,- p\®pI,. 

s,s' 

On the other hand, every state p^ ® p^, p^ € 'D{W), i = 1,2, is in 'D^{H^ (g) H^}, i.e. it is an extremal 
point of ViH^ (E) H^). Therefore it cannot be written as a non-trivial convex combination of elements 
from T>{Tl^ (g) H^), thus from a smaller set S{Ti^ (g H^). The invariance is obvious, since 

(Ti ® T2) o (pi P2) o {tI ® T2)t = (TipiT^) ® (Tip27^) 

and (TipiT^) e V\W) for pi e V^W). I 

Since we are working in a finite-dimensional space, the closeness of the corresponding hulls is automatic 
that can be derived from the following lemma. 

Lemma 5 IfVis an n- dimensional real vector space and x is a convex combination x = E^i ^i^i ^f 
certain points of V, then x is a convex combination of at most (n + l) points among Xi 's. 

Proof. It suffices to prove that x is a convex combination of (m — 1) of .t^'s, provided m > n + 1. Of 
course, we can assume that all ti > 0. If to > n + 1, then there are G R, not all equal 0, such that 
El* di = and E™ c-i^i = 0- There is io such that \aig/tig \ is maximal among \aig/tig\, i = 1, . . . , to. 
We can assume without loss of generality that io = to. Hence 

m-l , X 

a;= V t - ^^\x- 
and the above combination is convex, since iti — ^i^) > and 

m— 1 , m , m 

[ti - ) = 2_^[ti - J = } Ji = i. 
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Proposition 7 The convex hull conv{E) of a compact subset E in a finite dimensional real vector 
space V is compact. 

Proof. Suppose that the dimension of the space is n and denote by A„+i the compact (n + 1)- 
dimensional simplex 

n+l 

An+1 = {<= (ti,...,t„+i) : t, > 0,^t, = 1}. 

1 

According to the above lemma, conv(i?) is the image of the compact set A x x . . . x (i? appears 
in the product [n + l)-times) under the continuous map 

n+1 

A X E X . . . X E 3 {t,xi, . . .,Xn+i) i-> ^ t^Xi e V. 



Corollary 6 The set S{n^ ® H'^) is a compact subset ofu*{n^ H'^). 

The entangled states play an important role in quantum computing and one of main problems is to 
decide effectively whether a given composite state is entangled or not. An abstract measurement of 
entanglement can be based on the following observation (see also Ref. ,42') 

Let E be the set of all extreme points of a compact convex set K in a finite-dimensional real vector 
space V and let Eg be a compact subset of E with the convex hull Kq = conv(i?o) C K. For every 
non-negative function f : E ^ R+ define its extension Jk'-K^ R+ by 

fK{x)= inf VtJ(a.), (82) 

where the infimum is taken with respect to all expressions of x in the form of convex combinations 
of points from E. Recall that that, according to Krein-Milman theorem, K is the convex hull of its 
extreme points. 

Theorem 9 For every non-negative continuous Junction f : E ^ R+ which vanishes exactly on Eq 
the function fx is convex on K and vanishes exactly on Kq 

Proof. It is completely obvious that fx vanishes on the convex hull of Eq . The function fa is convex, 
since for every convex combination x ~ tiyi of points of K and every e > we can find extreme points 
Uj with convex combinations yi — s^aj and fxiyi) > ^Ifio^j) ~ £■ Hence 

fK{tiyi) = fK{tislaj) < t^sffiaj) < ti{f{yi) + e) = t^fxivi) + £■ 
Due to arbitrariness of e > we get 

fKitiVi) < t^fxiyi)- 

Note finally that fx vanishes exactly on Kq. Indeed Kq is compact due to propositionQand if a; ^ Kq, 
then x and Kq can be separated by a hyperplane, i.e. there is a linear functional </? : y R such 
that ^p{x) — a > Q and is negative on Kq. Denote by Ei the (compact) set of those points from 
E on which takes non-negative values and by F the minimum of f on Ei. Of course, -F > 0, since 
EiH Eq = 0. Let Af g R be the maximum of (p on E. Of course, M > 0. For any realization x = tiai 
of a; as a convex combination of points of E we have 

a = Lp{x) ='^ti(p{ai) < ^ tiip{ai) < M ^ t^. 

On the other hand, 

clF 

Y.^^f{(^^) > J2 *^/("') ^ ^ E ^ M ' 

so fxix) > ^ > 0. ■ 
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Corollary 7 Let F : V^{Ti,^ 'S'Ti.^) R+ be a continuous function which vanishes exactly onS^{Ti}® 
n^). Then 

M - ^i5CHi®w2) : V{n' ® H^) ^ R+ 

is a measure of entanglement, i.e. /i is convex and ^{x) = <^ x G Siji.^ (8) 7i^). Moreover, if f is 
taken U{Ti.^) x U{Ti.'^) -invariant, then is U{Ti.^) x U {Ti?)-invariant. 

Proof. The first part is a direct consequence of Theorem Also the invariance of ^ is clear: 

^l{TpT^) = inf(t,/(aO) = inf (i ./(Ta.Tt)) = inf(<,/(a,)) = 
w w w 

where T is in the corresponding group, 

W = {{t,, a^) : TpT^ = ^ a, G ^^(HI ® H^), > 0, ^ = 1}, 

and 



A careful study of the geometry of u*{{H^ ^Ti?) and criteria of entanglement we postpone to a separate 
paper. 
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